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Abstract. The strong-interaction shift ej* and broadening r™ a in pionic deuterium have been deter- 
mined in a high statistics study of the 7rD(3p — Is) X-ray transition using a high-resolution crys- 
tal spectrometer. The pionic deuterium shift will provide constraints for the pion-nucleon isospin 
scattering lengths extracted from measurements of shift and broadening in pionic hydrogen. The 
hadronic broadening is related to pion absorption and production at threshold. The results are 
e iP = (—2356 ± 31) meV (repulsive) and r^P = (1171^^) meV yielding for the complex 7rD scattering 
length a^ D = [-(24.99 ± 0.33) + i (6.22 + °2g)] x 10~ 3 m~\ From the imaginary part, the threshold 
parameter for pion production is obtained to be a — (251 \ t ) fib. This allows, in addition, and by using 
results from pion absorption in 3 He at threshold, the determination of the effective couplings go and gi for 
s-wave pion absorption on isoscalar and isovector NN pairs. 

PACS. 36.10.Gv, 25.80.Ls, 07.85.Nc Mesonic atoms, Pion inclusive scattering and absorption, X-ray spec- 
trometers 



1 Introduction 

Hadronic atoms reveal the influence of the strong force 
by a shift e and broadening r of the low-lying atomic 
levels with respect to the pure electromagnetic interaction. 
As atomic binding energies are negligibly small compared 
to the hadronic scale, a measurement of shift and width 
is equivalent to a scattering experiment at zero energy 
(threshold). Hence, such atomic data contain information 
on hadron- nucleus scattering lengths [TJ|2] . 

In pionic hydrogen, the atomic ground state level shift 
efj 1 and broadening r^ H are connected to the two isospin- 
separated pion-nucleon (ttN) scattering lengths. In the 
case of nuclei A(N, Z) with A > 2, eJ A and rf s A are 
measures of the complex pion-nucleus scattering length 
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a^A- The real part of o t a is attributed in leading order 
to scattering, whereas the imaginary part is due to pion 
absorption inducing itself a contribution to the shift [3J[3] . 

It has been shown, that strong-interaction effects in 
low Z pionic atoms can be related unambiguously to pion- 
nucleon (irN) scattering lengths as defined in modern field 
theoretical approaches of QCD [5] . Such threshold quan- 
tities are of great importance because they belong to the 
first scattering parameters accessible by lattice calcula- 
tions [6]. For the irN case, however, results are not yet 
available 0[S] . 

Pion-nucleus dynamics is understood to be built up 
from elementary 7rA — > ttN and irNN -H- AW processes 
taking into account nuclear structure, multiple scatter- 
ing, and absorptive phenomena. Vice versa, precise pion- 
nucleus data could set constraints on elementary ampli- 
tudes. Generally, systematic uncertainties involved in the 
extraction of a T A from atom data, if available, are smaller 
than in the extrapolation of cross-section data to thresh- 
old. In this case, normalisation problems cease to exist 
and corrections owing to Coulomb bound states are bet- 
ter under control than for scattering states (see sect. [3]). 
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Assuming isospin conservation and charge symmetry, 
ttN —> irN scattering at threshold is described completely 
by two real numbers, e.g., the isoscalar and isovector irN 
scattering lengths a + and a~ . They are defined in terms of 
two irN — > irN reactions or isospin 1=1/2 and 7 = 3/2 
by 

a+ = 2 ( a 7T-p^7r-p + a 7r+p^7r+p) = g ( a l/2 + 2a 3 / 2 ) , 

a =2 ( a TT-p^7T-p _ a 7T + p->7r + p) = g ( fl l/2 _ fl 3/2)- (l) 

Furthermore, among others the relations a^-n^r-K-n = 
a Tr+p^Tv+p and a^-p^^n = - \f2a~~ hold. 

The scattering lengths a + and a~ play a key role in 
the modern low-energy approach of QCD based on effec- 
tive field theories (EFT). A perturbative method — chiral 
perturbation theory (xPT) — has been developed, which 
is intimately related to the smallness of the pion's mass 
due to its underlying nature as a Goldstone boson [9] 
[TUIITTIIT2] . The interaction is calculated order by order 
as an expansion in the — compared to the hadronic scale 
of 1 GeV — "small" quantities momentum, pion mass, fine 
structure constant thus including the (isospin-breaking) 
electromagnetic interaction, and the light quark mass dif- 
ference uid — m u (strong isospin breaking) [13, 14, 15, 16, 
[T71HFIITO] . In such a theory short-range contributions are 
summarized in so called low-energy constants (LECs) which 
must be determined from ttN experiments (for details 
see [5J). 

From the scattering lengths a + and a - , important 
quantities of the ttN interaction are derived [5D] . The ttN 
a term jT2"ll2Tl[2"2"] is a measure of explicit chiral symme- 
try breaking, and the ttN coupling constant, obtained 
from the Goldberger-Miyazawa-Oehme sum rule [23 , 24 , 
[2"5] . is related to the strength of threshold photo pro- 
duction of charged pions [26 p7ll28] and the induced pseu- 
doscalar coupling g p as determined from muon capture in 
hydrogen [5511501151] . 

Access to a + and a~ is available from the Is level shift 
and broadening in pionic hydrogen (7rH): 

e? s H oc a n - p _^„- p = a + +a~~ + ... (2) 

IT? cc(l + ^j K- p ^o„) 2 cx (a-) 2 + ... . (3) 

Ellipses stand for corrections owing to the fact that the 
pionic atom constitutes a Coulomb state as well as for elec- 
tromagnetic and strong isospin and non-isospin breaking 
corrections, which are essential to extract a + and a~ in a 
well defined way from the data. These corrections may be 
calculated, e.g., within the framework of xPT and are of 
the order of a few per cent [51 [52"1I551154"1 I55] . 

In the 7rH case, the level width is as well due to scatter- 
ing — the charge exchange reaction n~p — > Tr°n — and, there- 
fore, is proportional to the square of a~ [53]. The Panof- 
sky ratio P = r(ir~p — s- n° n) / T '(tt~ p -> 771) = 1.546 ± 
0.009 37 accounts for the relative strength of radiative 
capture, the second channel contributing to the total broad- 
ening. For a recent analysis of 7rH atomic data see ref. [38] . 



It is known already from current algebra, that a + van- 
ishes in the chiral limit [9lll0j acquiring a finite value only 
because of chiral symmetry breaking [T2TI59"] . Hence, the 
s-wave ttN interaction is dominated by the isovector con- 
tribution with the consequence that it is difficult to extract 
a + from pion-nucleon and nucleus data with 7^0. 

Pionic deuterium (7rD) is of particular interest being 
the first natural choice to test the interaction of pions with 
isoscalar nuclear matter. Here, the leading order (LO) one- 
body contribution to the real part of the pion-deuteron 
scattering length a^d depends only on the isoscalar quan- 
tity a + . Considering the deuteron as a free proton and 
neutron plus corrections, one may write 

oc Rea^D 

2a+ + .... (4) 

Besides the corrections mentioned for the case of 7rH, el- 
lipses include here terms arising from multiple scattering, 
depending both on a + and a~, and absorptive contribu- 
tions. Noteworthy, that the second order correction, domi- 
nated by (a~) 2 , exceeds even the magnitude of the leading 
order term because of a + <C a~ [911T0]. Nuclear structure 
is taken into account by folding with the deuteron wave 
function [4"0"ll4Tll4"2"] ■ The exact relat ion between and 
Rea^D is discussed in detail in scct. |3.1| 

Hence, the isopin scattering lengths are accessible from 
the shift measurements in 7rH and in 7rD already with- 
out ultimately precise data for the hadronic broadening 
in hydrogen being available [23] . An effective constraint 
on a + and a~ , however, is obtained when combining the 
triple e^f, and e^f . A complete set of data has been 
achieved for the first time during the last decade [1511441 
E5JI46]. 

Such a constraint is highly desirable, because different 
methods exist to obtain the corrections needed for the ex- 
traction of the pure hadronic quantities a + and a~ from 
7rH and 7rD measurements. Recognising, among others, 
that corrections obtained using potential models were in- 
complete, significant effort went into calculations to estab- 
lish an unambiguous relation between atomic data and a + 
and a~ , in particular within the framework of xPT. It has 
been proven to be extendable in a defined way to 3-body 
interactions [17] and allows to treat electromagnetic and 
strong isopin-breaking terms on the same footing [5j|24, 48 , 
|4^1150115T115211551I5411551156115711581I591 . For a discussion and 
comparison of potential and EFT approaches see ref. [5J . 

As mentioned above, corrections for the ttN scattering 
lengths as obtained from 7rH amount to a few per cent 
of the leading contribution in complete contrast to the 
real part of the 7rD scattering length. Here, because of 
the smallness of a + , the isospin-breaking contribution was 
found to be as large as 40% in NLO xPT [58]. 

Finally, it turned out that for the extraction of the 
scattering lengths a + from the atomic data within the 
framework of xPT, the accuracy achievable is determined 
by one particular combination of LECs. It appears both in 
the correction term for e 7 ^ and , and thus constitutes 
a limit from the theoretical side [60] . In the case of 
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the LECs involved are known better [5,61 . Therefore, here 
the uncertainty of a~ is dominated by the experimental 
accuracy. 

Being different from the 7rH case, the hadronic broad- 
ening in 7rD and, consequently, the imaginary part of the 
a^D is not related to ttN scattering lengths a + and a~ but 
to the strength of s-wave pion absorption. Open channels 
in negative pion absorption in deuterium at threshold are 



n~d — > nn 
— > nnj 
— > nne + e~ 



— > nnn 



(5) 
(0) 
(7) 
(8) 



True absorption Q represents the inverse reaction of pion 
production in nucleon-nucleon collisions, i.e., Ima^D is 
also a measure of the strength of s-wave pion production at 
threshold. For that reason, pion-production cross sections 
were used to estimate the 7rD level width before it was 
experimentally accessible [5j!j • First attempts have been 
made to calculate the production strength rigorously also 
within the framework of xPT [6311 64] to be compared with 
the precise atomic data. 

The relative strength of the dominant channels true 
absorption and radiative capture ^ was measured to be 
S = r(ir-d -> nn)/r(TT~d -> mi) = 2.83 ± 0.04 [55]. 
The branching ratio of internal pair conversion ^ was 
determined to be 0.7% [66]. Charge exchange ^ is parity 
forbidden from s states and, therefore, only odd partial 
waves contribute which results in a fraction as small as 
(1.45±0.19) • 10 -4 [57]. The corresponding partial 3p level 
width can be neglected here. Hence, the relative strength 
of true absorption to all other final states contributing to 
the Is level width is found from the measured branching 
ratios to be 



S' 



r(nn) 



r(nnj) + r(nne + e ) 



2.76 ± 0.04, 



(9) 



i.e., about 2/3 of the hadronic width is related to the pion 
production/absorption process NN «-» irNN. The exact 
relation between pion-production strength and Ima^D is 
derived in sect. 13.21 

The aim of this experiment is to provide data on the 
7tD hadronic shift and width at least at about the level 
of accuracy achieved in recent or envisaged in ongoing 
theoretical calculations. A precise value for the shift will 
provide, together with the forthcoming new precision data 
for 7rH [68,69 , improved constraints for the irN scattering 
lengths a + and a - . As a first result of this experiment, the 
s-wave pion production strength derived from the width 
has been outlined recently [70] . 



2 Atomic cascade 

After pion capture in hydrogen isotopes, a quantum cas- 
cade starts from main quantum numbers n « 16 [71 , 72J 
(Fig.0. The upper and medium part of the atomic cas- 
cade is dominated by the collisional processes inducing 



Stark mixing, Coulomb de-excitation, and ionization of 
neighbouring atoms (external Auger effect). In the lower 
part, X-ray emission becomes more and more important. 
In light atoms, only electric-dipole transitions contribute 
preferring steps with maximal An. The cascade time, es- 
timated to be of the order of 0.1ns for densities around 
lObar |73j . is much shorter than the pion's life time. In ad- 
dition, molecular formation 7rD + D2 —¥ [(dd?:)d]ee must 
be considered. The influence of the above-mentioned pro- 
cesses is discussed in the following. 

Radiative de-excitation is the only process where col- 
lisional effects on the line energy and width can be ne- 
glected. Doppler broadening from thermal motion is of the 
order of 10 meV. Pressure broadening amounts at maxi- 
mum to 0.2 meV for the densities and temperatures used 
in this experiment. 

Stark mixing essentially determines the X-ray yields 
in exotic hydrogen. Because these systems are electrically 
neutral and small on the atomic scale, they penetrate sur- 
rounding atoms and, thus, experience a strong Coulomb 
field. Non-vanishing matrix elements (nlm\ E| nl'm!) in 
the presence of the electric field mix atomic states of the 
same principle quantum number n according to the selec- 
tion rules Al = ±1 and Am = [73]. An induced s state, 
from where pions disappear by nuclear reactions, leads to 
a depletion of the X-ray cascade. As the Stark mixing rate 
is proportional to the number of collisions during the ex- 
otic atom's life time, it explains the strong decrease of 
the X-ray yields with target density [73] . The yield of the 
7rD(3p — Is) X-ray transition has been measured at tar- 
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Fig. 1. De-excitation cascade in pionic deuterium. The en- 
ergy of the X-ray emitted in the 7rD(3p — Is) transitions is 
3.075 keV. Hadronic shift ei s and width 7\ s are of the order 
of 2.5 and 1 eV, respectively. The hadronic shift is defined as 
t\ s = E cxp — Eqed, i.e., a negative sign as is the case for 7rD 
corresponds to a repulsive interaction. 
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get densities equivalent to 15 bar and 40 bar STP to be 
(3.28 ± 0.43)% and (1.72 ± 0.16)%, respectively [75]. 

The natural line width of the atomic ground state tran- 
sition (3p — Is) of about leV is dominated by the life 
time of the Is state. Nuclear reactions from the 3p level 
contribute with less than 1 ueV and the 3p-level radiative 
width amounts to 28 fj,eV only. An estimate for the in- 
duced width from Stark mixing (3p O 3s) and external 
Auger effect at the 3p state, based on transition rates for 
7rH given in ref. [76], yields also negligibly small contribu- 
tions of < lueV at the target densities considered here. 
With a hadronic shift in 7rD being a factor of about 3 
smaller than in 7rH, Stark mixing might be enhanced, but 
strong effects are excluded because similar K yields have 
been observed for hydrogen and deuterium [75] . 

However, significant broadenings of the X-ray line in- 
duced by other cascade processes must be considered, in 
particular in the case of Coulomb de-excitation [75] • Here, 
the energy release during the transition (n — > n') is con- 
verted into kinetic energy of the collision partners being 
the excited 7rD system and, predominantly, one atom of 
an D 2 molecule [7611771175] . Such transitions dominate de- 
excitation from the start of the quantum cascade down to 
n « 10 and have been found to contribute even at lowest 
densities [50J . 

In lower-lying transitions, a significant energy gain oc- 
curs leading to a substantial Doppler broadening of sub- 
sequent X-ray transitions. The Doppler effect was directly 
observed in 7rH in time-of-flight spectra of monoenergetic 
neutrons from the charge exchange reaction at rest n~p — > 
7r°n 81, 55] and in an additional broadening of the line 
width of the muonic hydrogen (3p — Is) X-ray transi- 
tion 83; ■ The broadening turned out to be a superposition 
of several Doppler contributions attributed to different de- 
excitation steps n — > n' of the initial pTp system. Only 
An=l transitions could be identified. A similar behaviour 
was found for pionic hydrogen. Here, a significant increase 
of the total line width is observed with decreasing n of the 
initial state, which is attributed to the increasing energy 
gain of preceding Coulomb de-excitation steps [6T)1I54"] . 

Acceleration due to Coulomb de-excitation is counter- 
acted by elastic and inelastic scattering. This leads to a 
continuum of velocities below the ones well defined by 
a specific transition. Cascade calculations have been ex- 
tended to include the development of the velocity de- 
pendence during de-excitation and, therefore, are able to 
predict kinetic energy distributions at the time of X-ray 
emission from a specific level (extended standard cascade 
model ESCM [7611771175] ) . At present, calculations have 
been performed only for ttH [7^l771[7^[85l[86l[87] . Figure^ 
shows such an ESCM prediction for the 3p state. 

In the case of 7rD, the energies for the An = 1 Coulomb 
de-excitation transitions (7—6), (6 — 5), (5 — 4), and (4—3), 
are at 12, 20, 38, and 81 eV, respectively. Though kinetic 
energies differ only little from the 7rH case, velocities are 
significantly smaller (Fig.pJ possibly changing the balance 
between acceleration and slowing down. 

In the analysis of a measurement of the /iH(3p— Is) line 
shape it turned out [83], that first ESCM predictions are 



hardly able to describe the measured line width. There- 
fore, a model independent approach will be used to ex- 
tract the relative strength of Doppler contributions di- 
rectly from the X-ray data, which was applied success- 
fully first in the neutron time-of-flight [51] and also in the 
uH(3p — Is) analysis. Here, the continuous kinetic energy 
distribution is replaced by a few narrow intervals at the 
energies, where significant contributions are expected. Po- 
sitions and intensities of such contributions are then found 
by comparing a Monte-Carlo generated line shape to the 
data by means of a x 2 analysis. An approximation of the 
kinetic energy distribution by two narrow intervals, as 
used in the analysis of this experiment, is illustrated in 
fig. [2] A detailed description is given in sect. 5.2 



Auger transitions prefer de-excitation steps as small as 
possible, given that the energy gain exceeds the binding 
of the electron. Auger emission contributes mainly in the 
range An 6 — 10 with An = 1 transitions, but rapidly 
decreases for smaller n [7311771175] . Because of the small 
recoil energies, Auger emission cannot contribute to the 
high energetic part of the kinetic energy distribution. 

From muon-catalysed fusion experiments it is known 
that during /iH + H2 collisions metastable hybrid molecules 
are formed like [(ppfi)p]ee [551155] . An analogous process 
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Fig. 2. Prediction for the kinetic energy distribution of the 7rH 
atom at the instant of (3p— Is) transition for a density equiva- 
lent to 10 bar gas pressure together with the velocity scale. The 
velocities for the heavier 7rD system owing to the corresponding 
energy gains are given on the bottom. For the analysis of this 
7rD experiment, the energy distribution, after scaling the ener- 
gies from the 7rH to the 7rD case, was tentatively approximated 
by two narrow intervals representing a low-energy component 
(0-2 eV) and a second one for the high-energy component stem- 
ming from the (4 — 3) Coulomb de-excitation transition (for de- 
tails see sect. 5.2 1. The low-energy component dominates the 



distribution (note the broken vertical scale). 
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is expected in 7rD + D2 collisions [89] . Evidence that res- 
onant molecular formation takes also place from excited 
states has been found in muonic hydrogen [5U] . 

Such complex' are assumed to stabilise non radiatively 
by Auger emission. Possible X-ray transitions from weakly 
bound molecular states before stabilisation would falsify 
the value for the hadronic shift determined from the mea- 
sured X-ray energy due to satellite lines shifted to lower 
energies, whereas Auger stabilised molecules emit X-rays 
of at least 30 eV energy less than the undisturbed tran- 
sition |92| and are easily resolved in this experiment. Ra- 
diative decay out of such molecular states has been dis- 
cussed and its probability is predicted to increase with 
atomic mass [ST, 92 . A weak evidence for an Auger sta- 
bilised transition is discussed in sect. 16. II 

As molecular formation is collision induced, the frac- 
tion of formed complex molecules and the X-ray rate should 
depend on the target density. Hence, a measurement of 
the 7rD(3p — Is) X-ray energy was performed at differ- 
ent densities to identify such radiative contributions (see 
sect.Ol). 



3 Strong interaction 

In this section, the relations are introduced which will 
be used to extract t he c omplex pion-deuteron scatter- 
ing length a^D (scct. |6.5[ ) and t he t hreshold parameter 
for s-wave pion production (sect. 6.6 1 from the measured 



hadronic level shift and broadening. 



3.1 Scattering length 

The hadronic pion-nucleus scattering length a^A is related 
in leading order (LO) to the ns level shifts e ns and widths 
r ns by the Deser-Goldberger-Baumann-Thirring (DGBT) 
formula [TJ 



2aVc 4 

he 



Z 3 



,LO 
1 ttA 



(10) 



with fi being the reduced mass of the bound system of 
pion and nucleus A(Z,N). In eqs. (10 1 - (13) a denotes 
the fine structure constant. 

Shift and width are measured using a Coulomb bound 
state. Therefore, the complex quantity a^9 as extracted 
from eq. ( 10 1 is not identical with the pure hadronic scat- 



tering length a^A- Usually, the Coulomb interaction is 
taken into account by Trueman's formula P5}[M1[95] which 
expresses the correction by an expansion in the ratio scat- 
tering length-to-Bohr radius. In the case of hydrogen (Z = 
1), the complete expansion can be written up to order 
0(a 4 ) in a compact form yielding for the Is state of 

pionic deuterium |60j 



. r l6 



2a 3 /1 2 



he 

2a[ic 2 
he 



(11) 



(lna-l)-a 7rD + ^ ac 



The term <Sp ac =0.51% accounts for the interference of vacu- 
um polarisation and strong interaction |96j . Its uncertainty 
is assumed to be negligibly small compared to the exper- 
imental accuracy [5] . 

In the case of 7rD, the total correction to cl^d is only 
about 1% (see eqs. ([24]) and (25)). Hence, eq. (11 ) is solved 
with sufficient accuracy by insertin g th e leading order re- 
sult for a^g as obtained from eq. (10) and the influence 
of the experimental uncertainty must not be considered 
here. 

Real and imaginary part of the pure hadronic 7rD scat- 
tering length are then given by 



Re a^D — — 



he 



2a 3 /_tV 



(12) 



1 - 



Im cl^d 



he 

he r ls 



y mOL L ) Roa L 2 h0 D 



2aVc 4 2 



(13) 



(lna-l)-2Rea™ + 5™ c 



3.2 Threshold parameter a and pion production 

The cross section for the reaction pp — > ir + d is parametrised 
at low energies by [57] 



(14) 



where rj = p*^jm,^ is the reduced momentum of the pion 
in the wd rest frame. Here, a denotes the threshold pa- 
rameter representing pure s-wave pion production. Ap- 
proaching threshold, 77—5-0, higher partial waves (/3, ...) 
vanish. A value for a is obtained by fits to the quantity 
Up P ^. Tr+d/v an d then extrapolated to zero energy. In par- 
ticular at low energy, the correction factors Ci which take 
into account the Coulomb interaction are an important 
source of uncertainty in the determination. For example, 
the leading order correction Cg has been calculated to be 
as large as 30% with the additional difficulty to obtain an 
accurate error estimate |98U99j . 

A second experimental approach to a is to exploit the 
7rD ground state broadening, where uncertainties stem- 
ming from Coulomb correction factors and normalisation 
of cross sections are avoided. In order to derive in a model 
independent way the relation between a and Imo^D as 
obtained from pionic deuterium, purely hadronic (non- 
measurable) cross sections a are introduced to circumvent 
the problem of diverging Coulomb cross section at thresh- 
old. In this case Ci = 1, and a and (3 are pure hadronic 
quantities calculable, e.g., in the framework of xPT. The 
hadronic production cross section then reads 



r+d = ar) + Prj + 



(15) 



In addition, the reaction np —> ir°d can be used because in 
the limit of charge independence 2 • <r np ^ n a d 



'pp- 



Restricting to pion-deuteron s waves, in both processes the 
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same transition of a nucleon pair 3 Pi(I = 1) — » 3 S%(I = 
0) occurs, where true pion absorption nd nn induces 
the inverse reaction 3 Si[ 3 D 1 ](I = 0) -> 3 Pi(I = 1) on 
the deuteron's isospin nucleon-nucleon pair. 

Detailed balance relates pion production and absorp- 
tion by 



'VP- 



(16) 



with p* and p* being final state centre-of-mass (CMS) 
momenta |100| . Neglecting Coulomb and isospin breaking 
corrections, charge symmetry requires for the transitions 
7r~d — > nn and ir + d — > pp 



M. 



tt~ d—tnn 



1=1 M, 



TT+d—tpp 



(17) 



Isospin breaking effects are expected to be at most 1- 

2% [nmnD2| . 

A small difference in the transition rate results from 
the slightly larger phase space of the ir + d — >• pp reaction. 



&7T— d— >nn P 

d^rpp P 



V 



l \(s,ml,ml) 
A(s,m2,m2) 



0.982 . 



(18) 



The CMS momenta of proton and neutron, expressed in 
invariant variables using the triangle function A, are given 



by P* v 



J/2Vs [103]. The value corre- 



sponds to n ± d at threshold, where the total CMS energy 
squared is s = (m^ + m^±) 2 . The atomic binding energy 
of the 7r - D system is neglected. 

Combining optical theorem, charge invariance, detailed 
balance and inserting the parametrisation of the pp — > 
ir + d cross section (15), the imaginary part of the n~d — > 
nn scattering length reads in terms of the tt + production 
threshold parameter a 



Ima 



7T~d— >nn 



47T 
P% 
AlT 

_* 

P* 
4-7T 

1 

()7T 



1 7T~d— >nn 



&ir+d- 



2 

3 ' 

(p; 



p% 

Pn) 



T + d 



(19) 



To relate Imaj-^ given by pion production to the 
imaginary part of the pion-deuteron scattering length 
Ima n ]j as obtained from the atomic system 7rD, a cor- 
rection must be applied for the non true absorption chan- 
nels. Taking into account the relative strength 5" of true 
absorption to other s-state processes Q and exploiting 
A(x, y, y) — x(x — Ay) |103j . one obtains 



ImajD = (1 



= (1+1) 

= 2.48 • 10" 



Im a^-d-H 



v /(s-4m2)(s 



Ami, 



2Air m n 

l" 1 /Lt6 _1 • Q . 



(20) 



4 Experimental setup 

The experiment was performed at the 7rE5 channel of the 
proton accelerator at PSI, which provides a low-energy 
pion beam with intensities of up to a few 10 8 /s (Fig.p|. 
Pions of 112MeV/c were injected into the cyclotron trap 
II [10411105) and decelerated by using a set of degraders 
optimised to the number of pion stops in a cylindrical 
cryogenic target cell of 22 cm length and 5 cm in diameter 
placed in the centre of the trap. Its magnetic field, perpen- 
dicular to the beam direction, guides the particles after a 
few turns towards the target. 

The cell was filled with deuterium gas cooled by means 
of a cold finger to adjust the target density by tempera- 
ture. Data were taken at three different conditions, which 
are equivalent to gas densities of 3.3, 10, and 17.5 bar at a 
temperature of 20° C. The corresponding target parame- 
ters may be found in table[2j About 0.5% of the incoming 
pions are stopped per bar equivalent pressure of D2. X- 
radiation could exit the target cell axially through a 5 /^m 
thick mylar® window. The window foil is supported by 
1 mm thick horizontal aluminum bars constructed with a 
6 mm spacing. 

X-rays from the 7rD(3p— Is) transition were measured 
with a Johann-type Bragg spectrometer equipped with a 
spherically bent Si crystal cut along the (111) plane and 
having a large radius of curvature of R = 2982.2 ± 0.3 mm 
in order to minimise aberrations. The crystal plate was 
attached by molecular forces to a glass lens polished to 
optical quality, which defines the curvature. A possible 
miscut and its orientation were determined in a dedicated 
measurement [106 . The reflecting area of 10 cm in diam- 




ttE5 



Fig. 3. Setup of the 7rD experiment in the 7rE5 area at the 
Paul Scherrer Insitut (PSI). The roof of the concrete shield- 
ing is omitted to show the vacuum system connecting the cy- 
clotron trap (upper right), crystal chamber (upper left) and 
the cryostat of the X-ray detector (bottom left). Details of the 
mechanical setup are shown in figs.[4]and[5] 
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eter was restricted by a circular aperture to 95 mm in di- 
ameter to avoid edge effects and to 60 mm horizontally in 
order to keep the Johann broadening small [10711108] . 

Such a spectrometer is able to measure simultaneously 
an energy interval according to the width of the X-ray 
source when using a correspondingly extended X-ray de- 
tector. Being pixel detectors, charge-coupled devices 
(CCDs) are ideal detectors for X-rays in the few keV range 
because they combine an intrinsic position resolution with 
the good energy resolution of semiconductor detectors. In 
this setup an array of 3 x 2 CCDs was used covering in 
total 72 mm in height and 48 mm in width |109| . 

Monte Carlo studies show that about 2/3 of the inten- 
sity of the reflection is covered by the height of the CCD 
array. The solid angle of the crystal with respect to the 
target is 6 • 10~ 5 and the target fraction accepted by the 
crystal's angular width is s» 5- 10~ 3 . Including window ab- 
sorption and assuming a peak reflectivity of 40% for the 
Bragg crystal, the overall efficiency of the spectrometer 
results in about 10~ 7 . 

The CCDs are realised in open electrode architecture 
to minimise X-ray absorption in the surface structures and 
a depletion depth of about 30 ^m yields a quantum effi- 
ciency of about 80% in the 3-4 keV range. The pixel size of 
40 |um provides a sufficient two-dimensional position res- 
olution to measure precisely the shape of the diffraction 
image. The relative orientation of the individual CCDs 
as well as the pixel size at the operating temperature of 
— 100°C was obtained by means of an optical measure- 
ment using a nanometric grid [112] . The detector surface 
was oriented perpendicular to the direction crystal - de- 
tector. 

In the Johann setup, the energy calibration must be 
provided by a reference line of known energy. The best 



Table 1. Data to adjust the spectrometer in Johann setup. 
The Bragg angle Ob is calculated using a lattice distance d of 
2dm = (0.62712016± 0.00000001) nm [TTQ]. n denotes the or- 
der of reflection. For the 7rD(3p — Is) energy, here the experi- 
mental result of ref. [H] is given, yco = R ■ sin (Ob ± AOms) 
are calculated focal lengths taking into account the index of 
refraction shift A<9irs- The Ga Ka energies are taken from 
the compilation of ref. [Ill] . 





energy n Ob 


A0IRS 


ycD 




/eV 




/ mm 


Ga Kqi 


9251.674 ±0.066 3 39°52'22.3" 


2.4" 


1912.42 


?rD(3p - Is) 


3075.52 ±0.70 1 40°0'11.7" 


21.8" 


1917.84 


Ga K«2 


9224.484 ±0.027 3 40°0'44.1" 


2.4" 


1917.98 



angular matching for the 7rD(3p — Is) measurement was 
found with the gallium Ka fluorescence lines. Here, the Ga 
Kcv2 transition was chosen because of its smaller experi- 
mental error (Table[l}. The energy of the measurement 
line is obtained by the angular distance to the reference 
line, which is calculated from the position offset of the 
two reflections on the detector and the distance crystal - 
detector. 

The fluorescence target was made of a 25 x 20 mm 2 
GaAs plate mounted in the rear part of the gas cell out- 
side the pion stop volume and 82 mm off the centre of the 
cyclotron trap away from the crystal (Fig.|4|. X-rays were 
excited by means of an X-ray tube mounted on a window 
of the cyclotron trap chamber below the gas cell. 

In order to exclude the distortion of the 7rD(3]3 — Is) 
line shape from beam-induced excitation of Ga X-rays, 
data were taken using H2 gas at an equivalent density of 




Fig. 4. Sectional drawing of the connection of the cyclotron 
trap to the crystal chamber of the spectrometer. The symmetry 
axis of the magnetic field coincides with the line center of the 
trap to center of the Bragg crystal. 1: Bragg crystal, 2: traction 
relaxation, 3: concrete shielding, 4: gate valve, 5: X-ray tube 
(emission cone indicated), 6: position of pion beam without B 
field, 7: cryogenic target cell including calibration target (right 
end), 8: target cell support and cryogenic generator, 9: magnet 
coils of cyclotron trap, 10: LHe dewars. 



1 1 1 1 1 1 1 1 1 1 

Fig. 5. Sectional drawing of the connection of the crystal 
chamber to the X-ray detector. 1: concrete shielding, 2a and 
2b: inclination sensors, 3: Bragg crystal, 4: compensator, 5: lin- 
ear tables for crystal and detector adjustment, 6: gate valve, 
7: 5 (im Mylar window, 8: CCD array, 9: cold trap, 10: copper 
braid, 11: cold finger, 12: translation table for focal length ad- 
justment, 13: detector readout electronics, 14: turbomolecular 
pump, 15: air cushion support, 16: crystal chamber support. 
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10 bar. Using a hydrogen filling guarantees similar stop- 
ping conditions for the pion beam. No Ga X-rays could 
be identified during a measuring time corresponding to 
about 15% of the one used for D 2 at the same density. 
Other sources of background are discussed in sect. [5} 

The distance from the centre of the crystal to the 
centre of the cyclotron trap was 2100 mm, about 10% 
outside the Rowland circle given by the focal condition 
R ■ sinOft. The advantage of placing the X-ray source 
somewhat off the focal position is an averaging over non 
uniformities of the target. Both GaAs and D 2 target were 
extended enough that no cuts in the tails of the reflec- 
tion occur. The distance crystal - detector (Pig.j5J), cho- 
sen to be at the assumed 7rD focal length, was found to 
be ycD = 1918.1 ± 0.5 mm from a survey measurement. 

Alternating measurements of the Ga fluorescence ra- 
diation and the 7rD line were performed at least once per 
day. The Ga fluorescence X-rays were in addition used to 
monitor the stability of the line position. 

The spectrometer response was measured at the en- 
ergies 3104, 2765, and 2430 eV using the narrow Ml X- 
ray lines from helium-like argon, chlorine, and sulphur 
produced in a dedicated electron-cyclotron ion resonance 
trap (ECRIT) [IT^IlTiinTS] . It turned out that the res- 
olution function at a given energy can be built up from 
the ideal response calculated from the dynamical theory 
of diffraction for a perfect flat crystal (intrinsic resolu- 
tion) convoluted with the geometrical imaging at the mea- 
surement position by means of a Monte-Carlo ray-tracing 
code and by folding in an additional Gaussian contribu- 
tion. The intrinsic resolution is calculated here with the 
code XOP [116] . The Gaussian models possible imperfec- 
tions of the crystal material and mounting and was found 
to be sufficiently precise, that no difference is visible by 
eye between data and Monte-Carlo simulations. 

The value for the Gaussian width at the energy of the 
7rD(3p — Is) transition energy of 3075 eV was found from 
the fit to the results for argon, chlorine, and sulphur to 
be (122±8)meV. The narrow structure in fig. [8]- middle 
shows the Monte-Carlo generated response for the setup 
of the 7rD experiment. The total width of this resolution 
function is (436 ± 3) meV (FWHM) and close to the the- 
oretical limit of 403 meV — calculated by means of XOP — 
for the intrinsic resolution of a silicon crystal cut along 
the (111) plane. 

Details on the experimental setup and analysis may be 
found elsewhere [117] . 

5 Analysis 

Raw data of the X-ray detector consist of the digitized 
charge contents and a position index of the pixel. The 
granularity of the CCDs allows for efficient background 
rejection by means of pattern recognition (cluster anal- 
ysis). Photoelectrons from few keV X-ray conversion are 
stopped within a few micrometer only and, therefore, de- 
posit charge in one or two pixels with one common bound- 
ary. Beam induced background, mainly high energetic pho- 
tons from neutrons produced in pion absorption and cap- 



500 
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ADC/channels 



Fig. 6. Spectrum of the collected (digitized) charge in one 
CCD before and after cluster analysis of one of the 3 CCDs 
covered by the 7rD(3p— Is) reflection. Data are from the mea- 
surement at 10 bar equivalent density and correspond to 35% of 
the 7rD statistics collected with CCD 1. The peak around ADC 
channel 570 is due to the transition energy of 3.075 keV. Ver- 
tical dashed lines indicate the applied "energy cut" of ±2.5<r 
allowing for further background reduction. 



tured in surrounding nuclei, lead to larger structures. To- 
gether with a massive concrete shielding (Fig.[3| such events 
are highly suppressed. Defect pixels are masked by soft- 
ware. 

At first, the cluster analysis is performed. As expected 
only single (« 75%) or two pixel events (~ 25%) contribute. 
The spectra of the collected charge cleaned in this way 
show a pronounced peak originating from the 3 keV 7rD 
X-rays (Fig.|6|. For each CCD an individual energy cali- 
bration was performed by using the 7rD(3p— Is) line, nec- 
essary because of the different gain and noise behaviour 
of the various devices. The second calibration point is 
zero because the noise peak is suppressed on-line during 
data acquisition which reduces the amount of data to be 
recorded substantially. The energy resolution in terms of 
charge is determined by means of a Gaussian fit and found 
to be between 170 and 300 eV at 3keV. 

The good resolution in energy (i.e., in collected charge) 
by CCDs allows to apply a narrow "energy cut" in the 
ADC spectra (Fig. [6]). Thus an additional and significant 
background reduction is achieved in the two dimensional 
position spectra (Fig. [7]- top). Assuming a Gaussian shape 
for the energy resolution and setting various windows rang- 
ing from ± la to ± 4cr, the influence was studied of the 
peak-to-background ratio on the result for the hadronic 
broadening and found to be marginal. The minimum sta- 
tistical error is achieved for a cut of ±2.5er. The X-ray 
count rates at the equivalent densities of 3.3, 10, and 
17.5 bar are given in table[2J 

The hit pattern of the X-rays on the CCD surface 
shows a curvature originating from the imaging properties 
of the crystal spectrometer. The curvature is corrected by 
means of a parabola fit before projecting onto the axis of 
dispersion, which is equivalent to an energy axis (Fig. [7]). 
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The penetration depths of the 7tD(3p — Is) and the 
Ga X-rays differ significantly (5 and 105 fim for 3.075 
and 9.224keV, respectively [118J). Ga K X-rays also con- 
vert at or beyond the boundary of the depletion region 
of the CCD, where charge diffusion is already significant. 
Therefore, in this case compact clusters up to size 9 were 
accepted. As the Ga calibration measurements were per- 
formed without pion beam in the experimental area, no 
background events occur and, hence, no suppression algo- 
rithm is necessary to clean the Ga spectra. For one Ga 
calibration measurement, typically 4000 Ga Kct2 events 
were recorded. The total line width was found to be (2.9 ± 
0.1) eV, which is two standard deviations above the natu- 
ral line width of (2.66 ± 0.13) eV as given in ref. |119j . 





1800 




1600 




1400 




1200 


in 




OJ 

c 


1000 


c 




o 




u 


800 


\ 
> 






600 




400 




200 








)^fe:$ : i"7TD(3p-1s) 



800 850 

x/channels 



140- 






120- 




rrD(3p-1s) 


100- 


27K 1.1 bar 


Si (111) 


i „. 

o 

° 60- 


(10 bar) 


8 B =40' 


40- 






20- 






n. 







700 750 800 850 900 950 

x/channels 

Fig. 7. Top - Scatter plot of the 7rD(3p — Is) reflection mea- 
sured at a density equivalent to 10 bar (NTP). Its curvature 
is obtained by a parabola fit. Bottom - projection to the axis 
of dispersion after curvature correction. One channel (pixel) in 
the direction of dispersion (x) corresponds in first order reflec- 
tion to 76.402±0.001meV. 



5.1 Energy calibration 



The tabulated values for the Ga Ka energies [111) were ob- 
tained using the compound GaAs [120] - the same mate- 
rial used in this experiment, i.e., a possible chemical shift 
is irrelevant for this calibration. The position of the Ga 



Table 2. Experiment parameters. Rates (counts/Cb) are 
meant for a detector width of 14 mm and normalised to the 
integrated accelerator proton current (typically 2mA) for "en- 
ergy cut" of ± 2.5er. The equivalent density of the target gas is 
expressed as pressure value at a temperature of 20° C (NTP). 
The number of Ga calibration runs for each density is indicated 
in the last column. The backgound rate is density independent 
and scales linearly with the width of the "energy cut". The 
background given here stems from a detector area of 17.28 cm 2 
corresponding to 3 CCDs. 



D2 density target parameters 


counts 


no. of Ga 


(equivalent) T 


pressure 


7rD(3p- Is) 


calibrations 


/bar / K 


/ bar 


total / Cb 




3.3 ±0.1 33 


0.51 


1448±49 1.98±0.07 


10 


10.0 ±0.3 27 


1.09 


4010±74 6.40±0.12 


8 


17.5 ±0.4 25 


1.36 


4877±80 7.35±0.12 


7 


background 




1.64±0.11 





Ka2 line was determined applying a single Voigt profile in 
the fit, which is also the procedure used for the tabulated 



values [120) . The matching of the angular positions owing 
to the Ga Ka2 and 7rD(3p — Is) energies is evident from 

The 7rD line was modeled both with a Voigt profile, 
where Doppler broadening and response function together 
are approximated by a single Gaussian, and the true re- 
sponse as determined from the ECRIT data convoluted 
with the imaging properties by means of a ray tracing 
MC code and including the Doppler contributions from 
Coulomb de-excitation. Both methods yield the same po- 
sition value within a few hundreds of one CCD pixel. 

Various corrections to the measured line positions must 
be applied. Significant are, besides the index of refraction 
shift, bending and penetration correction, a shift of the 
gallium line due to the non-uniform illumination of the 
fluorescence target, and mechanical shifts from temper- 
ature changes from the nitrogen filling of the detector. 
These corrections are discussed in the following and listed 
together with less significant corrections in table[3] 

— Index of refraction. The Ga calibration line and the 
7rD(3p— Is) transition were measured in third and first 
order diffraction (table]!]). As the index of refraction 
strongly depends on the X-ray wave length a significant 
correction must be applied in order to obtain the right 
Bragg angle difference. The index of refraction A(9trs 
is calculated here from the atomic scattering factors 
as used by the code XOP [116] . The accuracy of such 
amplitudes is claimed to be of the order of 1% for 3 
and 9 keV, respectively, as are the angular corrections 
derived from them. The amplitudes of various sets vary 
by less than 0.5% [151111221 . 

— Crystal bending. Crystal bending leads to a depth 
dependence of the lattice spacing d, where close to the 
surface the maximum lattice distance is assumed. Be- 
cause of the difference in energy and accordingly the 
different mean penetration depth of the X-rays, the 
average lattice spacing and, consequently, the diffrac- 
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Fig. 8. Top - Ga Kot2 calibration line. The Kqi component 
is slightly reduced by the asymmetric emission characteristics 
of the X-ray tube. Middle - sum spectrum of the 7rD(3p — 
Is) transition measured at the equivalent densities of 10 and 
17.5 bar. The narrow structure represents the response function 
of the crystal spectrometer (normalised to the height of the 
measured spectrum). Bottom - residuum of data and fit. The 
fit is performed without assuming any Doppler components. 



tion angle is different. The change of d depends on the 
elastic properties of the nonisotropic crystal material 
and is calculated following the approach of ref. 123 
using the value v' = 0.182 for Poisson's ratio together 
with an expression for the angular correction as given 
by ref. [TH]. 

Penetration depth. The difference in penetration 
depth into the crystal leads to an offset of higher en- 
ergetic X-rays to smaller diffraction angles, which has 
been corrected according to the approach of ref. 11241 . 



The penetration depths (including absorption) were 
obtained from the code XOP to 0.72 //m and 3.84 /im 
for 3.075 and 9.225 keV, respectively. 



Fluorescence-target illumination. The intensity 
profile of the X-ray tube, used to illuminate the GaAs 
target, was found to decrease in the direction of disper- 
sion towards larger Bragg angles. Therefore, the centre 
of gravity of the line is shifted towards smaller Bragg 
angles or higher energies. The shape of the intensity 
distribution has been measured precisely and used for 
a Monte Carlo simulation to quantify the shift. 
Mechanical stability. Due to a fixed sequence for 
refilling the X-ray detector with liquid nitrogen (twice 
per day) and the gallium calibration (once per day af- 
ter filling) , small mechanical distortions arose from the 
cold nitrogen gas. These distortions were monitored 
every 15 min with two inclination sensors (of 1 /xrad 
precison), one attached to the crystal chamber and 
the second to the detector cryostat. In this way, both 
short and long term movements of the vertical crys- 
tal axis relative to the X-ray detector were recognized. 
For each of the three measurement periods an aver- 
age correction was determined from the inclinometer 
data. The maximal deviations were found to be about 
ztO.lmrad corresponding to about one pixel. 
Curvature correction. The curvature of the reflec- 
tion is determined from a parabola fit both to the 7rD 
and the Ga Kcv2 line, which can be assumed to be equal 
within the envisaged accuracy. The centres of grav- 
ity are determined choosing a left and right boundary 
aligned with the curvature of the reflection itself, the 
reflection divided vertically in 12 slices per CCD, and 
iterated until the minimal x 2 is reached. The maximum 
position by the various choices for the line and the 
boundaries are summarised in the given uncertainty. 
Defocussing. The small difference in the focal length 
of the 7rD and the Ga Ka2 leads to a small defo- 
cussing correction, which was determined by means of 
a Monte-Carlo simulation. 

Focal length. The parts constituting the mechanical 
connection between crystal and detector were adjusted 
to the assumed focal length of the 7rD(3p — Is) transi- 
tion. Their length and the non parallelity of the flanges 
were precisely surveyed resulting in an uncertainty of 
the distance of ±0.5 mm. 

Alignment. A possible deviation from the detector 
surface from being perpendicular to the direction crys- 
tal centre to detector centre leads to a correction for 
the measured distance of the 7rD and Ga Ka 2 reflec- 
tion. The uncertainty is composed of the setup of the 
CCD array inside the cryostat and the deviation from 
being parallel of the flanges of the connecting tubes, 
which were obtained from a survey measurement. 
Height alignment. If X-ray source, bragg crystal and 
detector centre are aligned perfectly in-plane, a sym- 
metric reflection with respect to this plane occurs at 
the detector surface. A misalignment causes a tilt of 
the reflection inducing a small error using the parabola 
approach for the curvature correction. The vertical po- 
sitions of target, crystal, and detector centre relative 
to the incoming pion beam were measured to be 206.5, 
205.1, and 206.1mm with an accuracy of ±0.2 mm. 
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— Pixel size. The average pixel size of the CCDs was 
determined in a separate experiment to be (39.9775 ± 
0.0006) /mi at the operating temperature of — 100°C 
[TT2] . 

Temperature normalisation. The crystal temper- 
ature was monitored regularly during the measure- 
ments. Though of minor impact, the lattice constant 
tabulated for 22°C, was readjusted to the average tem- 
perature of 30° C. The temperature was found to be 
stable within less than ± 1°C. Hence, any position vari- 
ation due to the change of the lattice constant can be 
neglected. 

As described above, the total uncertainty in the determi- 
nation of the 7rD(3p — Is) transition energy is dominated 
by the one of the Ga calibration line (± 27meV). The sys- 
tematic error of about 6m eV is about 50% of the statis- 
tical uncertainty (sect.pTJ) and comparable to the uncer- 
tainty of the calculated electromagnetic transition energy 
(sect.[6~2). 



5.2 Line width 

For comparison with data, the line shape is constructed 
by folding the response function (see sect. [4] and fig. 8]) 
with the Doppler induced width derived from a given ki- 
netic energy distribution (e.g., as shown in fig.^1), and a 
Lorentzian representing the natural width. Based on the 
experience of the study of the fiK(3p — Is) line shape [83"] . 
in a cascade model free approach, the kinetic energy dis- 
tribution is approximated by narrow intervals of a few 
eV width at energies which are inspired by the energy 
release of Coulomb de-excitation transitions. Such an en- 
ergy spectrum assuming only a low-energy (0-2 eV) and 
one high-energy component (owing to the (4 — 3) Coulomb 
de-excitation) is included in fig.[2j 

Spectra constructed in this way are compared to data 
by means of a x 2 analysis using the MINUIT package [125] . 
Natural line width, total intensity of the line, background, 
and relative weight of the Doppler contributions are free 
parameters of the fit. The sum of the Doppler contribu- 
tions is normalised to one. 

Following again the approach from the analysis of the 
fiH(5p — Is) transition, one tries to identify consecutively 
individual Doppler contributions by using first a single 
kinetic-energy component starting at energy zero and vari- 
able width. This corresponds to 7rD systems not being ac- 
celerated or already moderated down again by collisions 
to energies of a few eV. The \ 2 analysis shows that the 
upper boundary of the low-energy contribution must not 
exceed 8 eV. This result was achieved independently for 
the spectra taken at 10 bar and the 17.5 bar equivalent 
density (Fig. [9]). The result for the natural line width r ls 
turned out to be insensitive to the upper boundary when 
keeping its value at 8eV or below. Therefore, the low- 
energy component was fixed to the interval 0-2 eV in the 
further analysis. 

The most important higher energetic components are 
expected at about 80 and 115 eV stemming from the 



Table 3. Angular corrections (in seconds of arc (")) with as- 
sociated uncertainties as well as other sources of systematic 
and calibration errors occuring in the determination of the 
7rD(3p — Is) transition energy and hadronic shift ei s . The to- 
tal angular correction being the sum of the individual terms 
and associated uncertainties are converted to meV. The to- 
tal systematic error owing to experiment and setup (errors 
are assumed to be uncorrelated) is given conservatively as the 
quadratic sum and is calculated assuming on average the me- 
chanical stability of the 17.5 bar measurement. For comparison, 
the statistical error and the uncertainty of the QED transition 
energy (see sects. [tTT] and [rT2| is shown together with the one 
of the Ga calibration line. 



uncertainty 
correction 



7rD Ga Ka total correction 
/" /" /meV 



inrlp"v of rpfvapfinn 




2.4 


-42.6 


± 


0.4 


index of refraction 


21.8 




387.2 


± 


3.9 


crystal bending 


-3.73 


-3.65 


-1.4 


+ 


0.4 
0.9 


penetration depth 


-0.06 


-0.32 


4.6 


+ 


1.1 
2.7 


illumination GaAs target 





-1.075 


19.1 


± 


3.8 


mechanical stability 3.3 bar 


0.24 




4.2 


± 


2.5 


mechanical stability 10 bar 


0.71 




12.6 


± 


4.0 


mechanical stability 17.5 bar 


0.80 




14.2 


± 


1.3 


curvature correction 








± 


2.0 


defocussing 





-0.006 


0.11 


± 


0.01 


focal length 








± 


0.05 


alignment crystal-detector 








+ 




0.015 


height alignment 








+ 




0.016 


pixel size 








± 


0.14 


temperature renormalisation 








± 


0.02 


bias due to asymmetric errors 






-0.4 






total systematic error 








+ 


6.0 
6.6 


statistical error 








± 


10.9 


Eqed 








± 


7.9 


Ga Kq2 








± 


27 



Coulomb de-excitation transitions (4 — 3) and (5 — 3). 
Searches for any of these contributions failed, even when 
using the sum spectrum of the two measurements at 10 
and 17.5 bar. Varying energy or width of the high-energy 
Doppler contribution did not change this result. 

In this respect, pionic deuterium differs significantly 
from pionic and muonic hydrogen. There, such high-energy 
contributions were clearly identified from the width anal- 
ysis of the X-ray transitions [4l)ll6"9"l[53"] . For that reason, 
the ESCM calculation of the kinetic energy distribution 
for the 7rH(3p— Is) case, scaled to 7rD kinematics (Fig. [2]), 
is unable to reproduce the 7rD(3p— Is) line shape, because 
it contains rather strong contributions from (4 — 3) and 



12 



Th. Strauch et al: Pionic Deuterium 



499 



X 



XX> o O o 



□□ □ □ □ 



1083 



X 



5 10 15 20 

width low-energy component/eV 



25 



1080 



Fig. 9. Search for evidence of an extended low-energy com- 
ponent in the kinetic energy distribution by including in the 
fit of the data a uniform kinetic energy distribution of variable 
width starting at Tki n = 0. Allowing for kinetic energies above 
about 8eV increasingly downgrades the quality of the fit. Dia- 
monds (left x 2 scale) are due to the 10 bar data. Squares (right 
scale) are from a combined fit to the two spectra taken at 10 
and 17.5 bar. Degrees of freedom in the fits are 487 (left scale) 
and 974 (right). 



becomes more and more important with decreasing statis- 
tics. Where for the individual spectra measured at 3.3, 10, 
and 17.5 bar the bias was found to depend significantly on 
the statistics ((-38±3), (17±2), and (ll±2)meV), for the 
sum spectrum (10+17.5) bar it becomes almost negligible 
((-2±2) meV). The error of « 2 meV of the bias is given by 
the average of the results from the 400 Monte-Carlo spec- 
tra, which in turn determines the number of simulations 
necessary. 

The behaviour of the bias is understood as follows. 
For lower statistics, the suppression of the tails results 
in smaller widths, whereas in the case of higher statis- 
tics the fit tends to include background into the tail. The 
background has been found to be constant in all spectra. 
The systematical uncertainty of /i s due to different back- 
ground levels within the limits as obtained from the fit to 
the data is estimated to below 1 meV by means of Monte 
Carlo simulations. 

Again 400 Monte-Carlo generated spectra, based on a 
calculated ESCM energy distribution, were used to quan- 
tify the sensitivity to any further high-energy component. 



(5 — 3) Coulomb de-excitation. At present, no explanation 
has been found for such an unequal behaviour. 

It has been studied in detail which fraction of high- 
energy components may be missed in the fit by studying 
Monte-Carlo generated spectra for the statistics collected 
for the sum of the 10 and 17.5 bar measurements. The in- 
tensity of the Doppler contributions was not restricted to 
positive value s to allow an unbiased search for the mini- 
mum x 2 (Fig- 10 ~~ top). The normalization to one is then 



maintained by a correspondingly increased value for the 
other component. 

The probability to miss a Doppler component, stem- 
ming from kinetic energies around 80 eV and correspond- 
ing to the (4—3) transition, is displayed in fig. 10 - bottom 



as a function of its relative strength. It can be seen, that 
a contribution of 25% or larger can hardly be missed. For 
10% relative intensity, the chance is about 15% to find in- 
tensities < in the fit. Taking symmetric limits a rou nd the 
maximum at about 10% relative frequency (Fig. 10 - top) 



corresponds to la with respect to the full distribution. For 
each set of conditions 400 simulations were performed. 

Assuming no Doppler components yields an upper limit 
for ri s which is identical to the result obtained when using 
only the low-energy component of 0-2 eV. The residuum 
of such a fit is shown in fig. [8]- bottom. The limit of sen- 
sitivity for a component at 80 eV kinetic energy of 10% 
yields a lower bound for J~i s (— A r sys ) corresponding to 



the above-mentioned la criterium (Fig. 10). The distribu- 



tion of the results for the weight of this component and 
the Lorentz width r reflects the fluctuations owing to the 
limited statistics (Fig. 111. 



From the various above-mentioned sets, each with 400 
simulations, a possible systematic deviation (bias) of the 
analysis code was examined. Such a deviation stems from 
an imperfect description of the probability distribution 
by the fit model 126 , here used to extract ri s , which 




25 



relative intensity high-energy component 



MC 



— <>- 



0.05 0.10 0.15 0.20 0.25 0.30 

relative intensity high-energy component 

Fig. 10. Top - distributions of relative intensity of a Doppler 
broadened component induced by the (4 — 3) Coulomb de- 
excitation transition found from fits to Monte-Carlo (MC) gen- 
erated spectra. The relative strength used for the simulations is 
10% (filled) and 25% (hatched) relative intensity, respectively. 
In both cases, the result is based on 400 spectra generated for 
the statistics of the sum of the 10 and 17.5 bar data. Bottom 
- probability to miss a Doppler contribution corresponding to 
the (4 — 3) Coulomb de-excitation (modeled by a uniform ki- 
netic energy distribution from (76 — 84) eV) as a function of its 
relative intensity. 
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An additional 6% fraction corresponding to the predicted 
strength of the (5 — 3) de-excitation together with a 10% 
contribution from the (4 — 3) is identified by the fit in 
more than 2/3 of all cases. In the other 1/3 of the fits, the 
(5 — 3) fraction is absorbed into the (4 — 3) contribution 
leading to an overestimate of its intensity and a shift to 
higher energies. 

1400 1 1 

MC input 10% 

1300 

"A Tsys 

> 1200 . ■ . 

e . ^^^^is-:.r . 

1 1 oo 

1000 i i - .— -T-" 

900 J ■ 1 ■ 1 ■ ■ ■ ■ 

-0.1 0.1 0.2 0.3 

relative intensity high-energy component 

Fig. 11. Distributions of the Lorentz width r versus the rela- 
tive intensity of the Doppler broadened component stemming 
from the (4 — 3) Coulomb de-excitation. The results are ob- 
tained from 400 Monte-Carlo simulations performed for the 
statistics of the sum of the 10 and 17.5 bar data. Input values 
are 10% relative intensity for the Doppler affected component 
and 1130eV for T. 



A weighted average is calculated for the 7rD(3p — Is) 
transition energy from the results for the three different 
target densities by using the individual contributions to 
the statistical and systematical error. The final value and 
error (quadratic sum) is obtained when combining these 
errors with the uncertainty of AEc a for the Ga Ka2 line: 

Kr>(3 P -u) = (3075.583 ± 0.030) eV. (21) 

The experimental uncertainty is dominated by the error of 
the Ga Ka 2 calibration line. In case a significant theoreti- 
cal progress requires further improvement, the calibration 
error could be reduced by a factor of about 2 in a dedicated 
measurement using a double flat crystal spectrometer. 



Table 4. Measured 7rD(3p — Is) transition energy and associ- 
ated uncertainties. 



equivalent 
density 
/bar 


E(3p - Is) 
/eV 


AE(3p - Is) 
stat sys AEgh 
/ meV / meV / meV 


3.3 


3075.509 


±28 


+ s ±27 


10 


3075.594 


±17 


+_l ±27 


17.5 


3075.599 


±16 


+ J; ±27 


weighted average 


3075.583 


±11 


±7 ±27 



No evidence for an additional, i.e., asymmetric broad- 
ening from non resolved satellites due to molecular for- 
mation has been found within the experimental accuracy. 
This is corroborated by the observation that a possible 
density dependence of the 7rD(3p — Is) en ergy contradicts 
the expected behaviour (see sects. 6.1 and 6.3 1. 



6 Results and discussion 

6.1 Experimental 7rD(3p — Is) transition energy 

The results for the 7rD(3p — Is) X-ray energy are summa- 
rized in table[4]taking already into account the corrections 
and their uncertainties as given in table[3j The values ob- 
tained for the three different target densities are consistent 
within two standard deviations (neglecting the common 
uncertainty of the Ga Ka 2 energy). 

It is concluded that we cannot identify radiative de- 
excitation from molecular states within the experimental 
accuracy, even more, as it is expected that such a process 
decreases the 7rD transition energy with density. Notewor- 
thy, that also the analysis of the /iH(3p — Is) line shape 
yields no evidence for a broadening owing to molecular 
effects [83] . 

Inspecting the background at the low-energy side of 
the 7rD transition for possible Auger stabilised molecules, 
a weak structure was found separated by (28.6 ± 0.5) eV. 
Its relative intensity was determined to be (3.4 ± 1.6) -10 -3 
of the main line. 



6.2 Electromagnetic 7rD(3p Is) transition energy 

In order to extract the hadronic shift, the pure electro- 
magnetic transition energy has been recalculated within 
the framework of quantum electrodynamics (QED). The 
contributions to the level energies are given in table [5] 
The error of the calculated QED transition energy Eqed 
= 3077.939 ± 0.008 eV is dominated by the uncertainty 
of the charged pion mass |127j (±7.72meV), whereas the 
ones of the charge radii of the deuteron (± 1.19 meV) and 
the pion [127 j (± 0.95 meV) contribute only marginal. The 
numerical accuracy is better than 1 meV. For this calcula- 
tion, a deuteron charge radius rrj = (2.12809 ± 0.00031) fm 
has been used as derived from the most recent value of 



the proton charge radius |128j and the radii difference ob- 
tained from the hydrogen-deuterium isotope shift mea- 
sured by means of two-photon spectroscopy [129] . 

The correction to the ground-state energy arising from 
the electric polarizability of the deuteron of (-26.3±0.5) 
meV is obtained by scaling from a calculation for muonic 



deuterium |130|I131] . Contributions from the polarizability 
of the pion |132j are three orders of magnitude smaller and, 
consequently, are neglected. 

The strong-interaction shift of the 7rD 3p state is esti- 
mated as for the s-wave in leading order by the sum of the 
pion-proton and pion-neutron interaction. The ird scatter- 
ing volume may be expressed by twice the isoscalar and 
angular momentum averaged parameter cq of the Kisslingcr 
potential 3,48,62 , which can be calculated from the irN 



14 



Th. Strauch et al: Pionic Deuterium 



Table 5. Contributions to the pure electromagnetic binding 
and the (3p— Is) transition energies. The terms vac. pol. 11, 13, 
and 21 stand for Uehling, Wichmann-Kroll, and Kallen-Sabry 
contributions. The shift value due to the polarizability of the 
deuteron is derived from results given in refs. [13011131] . A "0" 
indicates a negligibly small but finite value. 



-2.2 





3p 


Is 


E QED (3p- Is) 




/eV 


/eV 


/eV 


Coulomb 


-384.31079 


-3458.52422 


3074.21344 


self-energy + finite size 





0.00220 


-0.00220 


vac. pol. 11 


-0.01432 


-3.72948 


3.71516 


muon vac. pol. 11 





-0.00034 


-0.00034 


vac. pol. 13 





0.00002 


-0.00002 


vac. pol. 21 


-0.00013 


-0.02793 


0.02780 


recoil 1 


-0.00004 


-0.00297 


0.00293 


relativistic recoil 





0.00029 


-0.00029 


two-loop vac. pol. 





-0.00577 


0.00577 


pion finite size 





-0.04764 


0.04764 


deuteron polarizability 





-0.0263 


0.0263 




-0.000008 




0.000008 


atom recoil 






-0.00235 


total 


-384.32528 


-3462.26684 


3077.93922 



scattering volumes. Values are given, e.g., in ref. [133J. In- 
serting the result for 2 • cq in the Trueman formula for 
the 2>p state [9"3llM] , one obtains a negligibly small shift of 
e 3p ps 8 fieV. 



6.3 Strong-interaction shift e± s 

Subtracting measured and calculated QED transition en- 
ergies, the hadronic shift is obtained to 



(2356 ± 31) meV. 



(22) 



The negative sign indicates a repulsive interaction, i.e., 
the strong interaction reduces the atomic Is level binding 
energy. 

Our result is compared to previous measurements in 
table[7] and fig. [12} Note that the electromagnetic transi- 
tion energy Eqed differs slightly in the analysis of the 
various experiments. Worth mentioning, that the results 
of the three precision experiments are in good agreement 
though three different calibration lines have been used. 

Combining the results of this experiment and the ones 
given in refs. [H] and |3S] yields a weighted average of 
ei s — (— 2379± 19) meV assuming any absence of a density 
dependence for the 7rD(3p— Is) transition energy. 

Inspecting fig.[l"2"|suggests that an extrapolation to den- 
sity zero yields the value aimed at. A linear fit yields a 
slope of (6.3 ± 1.7) meV/bar and e u = (2445 ±27) meV at 
density zero. However, at present there is no explanation 
for such a density dependence (see sect. [2]), which might be 
subject for a next generation experiment and, therefore, 
we stick to the density independent average given in (22). 
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Fig. 12. Experimental information on the hadronic shift in 
pionic deuterium. The error of the individual measurements 
at 3.3, 10, and 17.5 bar equivalent density and, therefore, also 
of the weighted average is dominated by the uncertainty of 
the Ga K«2 energy. The results from the two previous experi- 
ments at 2.5 bar [35] and 15 bar [33] are adjusted according to 
the new values for the pure electromagnetic transition energies 
(Tabled . 



6.4 Strong-interaction width ri s 

Corresponding to the discussion in sect. |5.2[ the hadronic 
broadening is extracted by introducing only a low-energy 
component from 0— 2 eV in the kinetic energy distribution. 
The combined result is obtained by averaging according to 
the statistical weig ht (Table[6]). The large negative error is 
dominated by the uncertainty that a Doppler contribution 
of the level of about 10% might be missed in the analysis. 
Combining all errors (quadratic sum), wc obtain 



+ 23 

r u = ( 1171 _ 4g | nu\. 



(23) 



The result is in good agreement with the earlier mea- 
surements, but a factor of about 3 more precise (Tabled 



and fig. 13|). The weighted average of the results of this 

is r u = (i 



experiment, 



and 



1165 



+ 22 

- 38 



meV. 



Table 6. Pionic deuterium ground state width _Ti s as obtained 
from the individual measurements at 3.3, 10, and 17.5 bar and 
the sum spectrum at (10+17.5) bar. The weighted average is 
calculated from the result of the sum spectrum (10+17.5) bar 
and the one of the 3.3 bar measurement. 



equivalent 
density 
/ bar 


A. 
/meV 


Ar ls 

stat 
/meV 


sys 
/ meV 


3.3 


1246 


±71 


+ 
- 152 


10 


1177 


±38 


+ 
- 92 


17.5 


1121 


±32 


+ 
- 81 


(10+17.5) 


1162 


±24 


+ 
- 43 


weighted average 


1171 


±23 


+ 
- 43 
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Table 7. Transition energies and hadronic effects found for pionic deuterium and corresponding experimental conditions. The 
total error of the measured transition energy E exp is calculated quadratically from the statistical and systematical contributions. 
The hadronic shift is defined by ei s = E cxp — Eqed- Note the change of Eqed (*) due to a new calculation performed for the 
analysis of this experiment which shifts the result of [44 to ei s = —2419 ± 100 meV. Also for the (2p — Is) transition a new 
QED value of 2597. 519±0. 008 eV is obtained resulting in ei s = —2461 ± 55meV for the experiment of [35]. No uncertainties for 
the QED values have been reported for the experiments decribed in refs. 45 , 135, 136 . 



transition equivalent energy E exp Eqed eis As 

density calibration 

/bar /eV [eV /mcV /meV 



nD{2p- 


15) 


4-8.5 


Bi My edge 


2592.8 


+ 


1.6 
2.0 


2597.61 ± 0.15 


-4800 


+ 


1600 
2000 










nD(2p- 


Is) 


» 3.5 


Cu Ka 


2593.3 


± 


2.3 


2598.1 


-4800 


± 


2300 








135 


■wD(2p- 


Is) 


« 50 


CI Ka 


2592.1 


± 


0.9 


2597.6 


-5500 


± 


900 








[136] 


irD(3p- 


Is) 


15 


Ar Ka 


3075.52 


± 


0.07 


3077.95 ± 0.01* 


-2430 


± 


100 


1020 


± 


210 


EH 


irD{2p- 


Is) 


2.5 


CI Ka 


2595.058 


± 


0.055 


2597.527* 


-2469 


± 


55 


1194 


± 


105 


m 


7rD(3p- 


Is) 


3.3/10/17.5 


Ga K«2 


3075.583 


± 


0.030 


3077.939 ± 0.008* 


-2356 


± 


31 


1171 


+ 


23 
49 


this exp. 



Again, a possible density dependence might be sus- tering length reads 
pected (Fig|l~3"|). Performing a linear fit — taking into ac- 
count the asymmetry of the errors — yields As = (1229 ± 
75) meV at density zero and a slope of (— 6.8±4.8) meV/bar 
Also here, no reason was found for such a density depen- 
dence and we keep the density independent result. 
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Fig. 13. Experimental information on the hadronic broadening 
in pionic deuterium. Previous measurements were performed at 
2.5 bar [45] and 15 bar [33] equivalent density. 



6.5 Scattering length a^o 



Real and imaginary part of the ttD scattering length a ff rj 
as determined from the various measurements and includ- 
ing the Trueman correction are compared in table|8] 

The 7rD scattering length was studied recently within 
the xPT approach [60] . In this work, almost prefect over- 
lap is achieved of all constraints imposed on the irN iso- 
scalar and isovector scattering length a + and a~ by pio- 
nic hydrogen and deuterium data. In addition, in a com- 
bined analysis using the preliminary results for shift and 
width in 7rH as given in |69j together with the new shift 
value in irD from this experiment, the uncertainty for a + 
could be reduced by a factor of two yielding a positive 
value by two standard deviations. By using the results 
for a + and a~ from this combined analysis, Rea^rj = 
(-25.4 ± 2.6) x 10~ 3 m" 1 is obtained. 

The imaginary part was calculated in LO xPT to be 
Ima^o = (5.65 ± 1.75) ■ 10~ 3 to" 1 [53]. One result of this 
work is that the real part induced by this absorptive con- 
tributions turns out to be (7±2)% of the total Ima^D- It is 



much smaller than the usual naive estimate of Re a ™& uced 



Im a n a [1] , but can be traced back in the case of 7rD 
cancellations of individually sizeable terms. 



to 



Evaluating the DGBT formula (j 1 0|) , one obtains for the 
:at 

0.26)1- lO" 3 to; 



complex scattering length a^g = [[^25.07+0.33) +i (6.23± 



eqs. (12) and (131 one obtains 



Inserting these leading order values in 



Rea^D = 0.010642 m~ l /eV ■ 0.9968 
Ima^o = 0.010642 m^/eV- 0.9989 



(As/2) 



(24) 
(25) 



with the factors 0.9968 and 0.9989 representing the "True- 
man" correction. Consequently, our result for the 7rD scat- 



6.6 Threshold parameter a in pion absorption 

The value derived from the ttD hadronic width As reads 



a = 251 



-11 



fib . 



(27) 



The central value differs by 1 fib from the one given in 
ref. [7D] , because there the correction term 5™ c (see sect 
had not been considered. 
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Table 8. Real and imaginary part Rea^D and Ima-n-D of the 7rD scattering length as extracted from atomic data including the 
Trueman correction (12 1 and (13 1. Threshold parameter a for pion production and lma,D are corrected for radiative capture 
and pair conversion ( 20[ ) but omitting any isospin breaking contributions. Production experiments reporting only the statistical 
uncertainty are marked (*). Arrows (<— / —>) indicate primary and derived quantities. The last line shows the result of a 
combined analysis of 7rH and ttD data (see sect.|6.5[). Theoretical approaches are listed in chronological order. 
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Faddeev (Bryan-Scott) 
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Coupled channels 
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Heavy meson exchange 












203 


± 


21 


[14911150] 


Relativistic field theory 


- 57 
















nn] 


X PT NLO 






5.46 


± 


1.74 «- 


- 220 


± 


70 


[63] 


X PT NLO 






5.65 


± 


1.60 - 


-> 228 


± 


65 


64 



xPT analysis combining 7rH and 7rD 



25.4 ± 2.6 



60 



The theoretical understanding of NN f-> NNir reac- 
tion is continuously increasing. Within the approach of 
xPT, a study of pion production including the terms in 
NLO yields cv NLO = (220 ± 70) fib [53] in good agreement 
with the pionic deuterium results (Tableland fig 14). The 
theoretical uncertainty is expected to decrease to below 
10% within a few years from forthcoming NNLO calcula- 
tions [137] . 



The parameter a when determined from pion produc- 
tion experiments shows wide fluctuations even when com- 
paring recent data. Often, only the statistical error is given 
for the cross section of the production experiments, but 
the fluctuations suggest significant systematic uncertain- 
ties, which may arise from uncertainties in the normali- 
sation and/or Coulomb corrections. The expected order 
of magnitude of isospin breaking effects is about or less 
the precision of this experiment, which is far below the 
variation of the pion production data. 



6.7 Pion absorption on nucleon-nucleon pairs 

To describe the s-wave absorption strength on isoscalar 
(1=0) and isovector (1=1) nucleon-nucleon pairs (NN), 
effective couplings go and g\ have been introduced rep- 
resenting the transitions 3 Si(I = 0) — > 3 P\(I = 1) (go) 
and 1 So(I = 1) 3 Po(I = 1) (ffi) BSD- The transition 
strength | go 1 2 is hence related to the imaginary part of 
the scattering length a T D as determined from Z~i s in nD by 
ImciTrD = | (p'n) 3 I go | 2 (1+1/S") [153] and, consequently, 
also to the threshold parameter a by 



30 | 2 = 



1 



27rm, 



Pp 

P'n P 



1 



(28) 



(see sect 



3.2|. The CMS momentum of the neutrons ac- 

2.6076 m„. 



quired in the absorption reaction is p* 

The relative strength of the coupling for absorption on 
isoscalar to isovector NN pairs was obtained to | 50/31 |= 
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a previous pionic-deuterium data 
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v theoretical approaches 



Fig. 14. Threshold parameter a for s-wave pion production 
compared to results from previous experiments and theoreti- 
cal studies(from left). Previous X-ray measurements: [45], |44| : 
pion production and absorption: [T3S], [135], [HD], [HI], EH], 
[143] ,[144], [145]; Theoretical calculations: [97], [146], [98], 
[T47] , [I47j . [US] , [T4W5D] . [S3], [53]. Points are given in the 
same order as listed in tabled] 



2.48 ± 0.24 from the ratio of neutron-neutron to neutron- 
protons pairs emitted back-to-back after pion absorption 
in the helium isotopes 154, 155 . From a as given by Im a X D 
and \ go/9i I one obtains 



90 I =(2.780 nnes ]-10-"m 



(29) 



91 



+ 0.027\ 
- 0.058 , 

(1.12 ± 0.11) • 10~ 2 m- 2 , (30) 



7 Summary 

The 7rD(3p— Is) X-ray transition in pionic deuterium has 
been studied using a high-resolution crystal spectrometer 
to determine the ground-state strong-interaction effects. 

The accuracy of 1.3%, obtained for the shift, super- 
sedes the theoretical uncertainty of 9% achieved recently, 
and which can be traced back mainly to insufficient know- 
ledge of low-energy constants. However, when used as a 
constraint for the pion-nucleon scattering lengths a + and 
a~ as determined from pionic hydrogen, the new 7rD re- 
sult reduces the uncertainty for the isoscalar scattering 
length a + by a factor of about two yielding a non-zero 
and positive value [60] . 

The accuracy achieved for the hadronic broadening of 
4.2% reaches already the expected final uncertainty of (5— 
10)% of forthcoming NNLO %PT calculations [137 . From 
its value, the transition strength a for s-wave pion pro- 
duction in nucleon-nucleon collisions is determined with 
unprecedented accuracy. 

No line broadening by radiative de-excitation from non 
stabilised molecular states was identified which is corrob- 
orated from the study of the energy dependence of the 
7rD(3p — Is) transition energy. However, weak evidence 
(2<t ) was found for a Auger stabilised state. 

Noteworthy, that at the 10% level no high-energy com- 
ponents could be identified stemming from low-lying 
Coulomb de-excitation transitions in contrast to pionic 
and muonic hydrogen. 
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where the effective coupling g\ describes, e.g., the thresh- 
old transition strength of the reaction pp —> ppir [156) . 

A phenomenological analysis for NN — > NNir close 
to threshold using the ansatz a = Bog 2 |97] yields from 
fits to the cross section Bo = 27.6 /ib. The uncertainty is 
estimated to be as large as 50% |157j . However, the simple 
energy dependence is only valid in the limit of an infinitely 
large scattering length. 

Scaling our result for g\ one obtains a smaller value 
of Bo = 13.3 ± 1.4 /ib. This may be due to the fact that 
the extraction on both | 30/51 | from ir~He absorption 
and go from 7rD are affected differently by NN final state 
interactions. 

As discussed, e.g. bv |158([159j . the phenomenological 
ansatz neglects the strong NN final-state interaction lead- 
ing to a completely different and complicated energy de- 
pendence close to threshold, in particular, when compar- 
ing pn and pp ( 1 5'o) final states [160U161J . This prevents a 
direct determination of the production amplitude from a 
simple parametrization of the cross section like cx g 2 . 

Calculations including final-state and Coulomb inter- 
action lead to values for a/g 2 as small as 5/ib when ap- 
proaching threshold |158|IT59"] . This is in agreement with 
the long known suppression of the pp — > ppir amplitude 
in LO [146] and a %PT calculation showing that loop con- 
tributions vanish up to NLO |162j . 



We would like to thank N. Dolfus, B. Leoni, L. Stohwasser, 
and K.-P. Wieder for the technical assistance and V. Baru, 
C. Hanhart, A. Nogga, M. Hoferichter, and A. Rusetsky for 
continuous exchange of information on theoretical progress 
on the 7rD system. The Bragg crystal was manufactured 
by Carl Zeiss AG, Oberkochen, Germany. Partial fund- 
ing and travel support was granted by FCT (Lisbon) and 
FEDER (Grant No. SFRH/BD/18979/2004 and project 
PTDC/FIS/102110/2008) and the Germaine de Stael ex- 
change program. LKB is Unite Mixte de Recherche du 
CNRS, de l'Ecole Normale Superieure et de UPMC No. 
8552. This work is part of the PhD thesis of one of us 
(Th. S., Univ. of Cologne, 2009). 



References 

1. S. Deser, M. L. Goldberger, K. Baumann, and 
W. Thirring, Phys. Rev. 96, 774 (1954). 

2. D. Gotta, Prog. Part. Nucl. Phys. 52, 133 (2004). 

3. T. E. O. Ericson and W. Weise, Pions and Nuclei, (Claren- 
don, Oxford, 1988) chapter 6. 

4. K. A. Bruckner, Phys. Rev. 98, 769 (1955). 

5. J. Gasser, V. E. Lyubovitskij, and A. Rusetsky, Phys. Rep. 
456, 167 (2008). 

6. M.Liischer, Commun. Math. Phys 105, 153 (1986); Nucl. 
Phys. B 354, 531 (1991). 



18 



Th. Strauch et al: Pionic Deuterium 



7. A.Torok et al, Phys. Rev. D 81, 074506 (2010). 

8. S. R. Beane, W. Detmold, K. Originos, and M. J. Savage, 
Prog. Part. Nucl. Phys. 66, 1 (2011). 

9. S.Weinberg, Phys. Rev. Lett. 17, 616 (1966). 

10. Y.Tomozawa, Nuovo Cim. A 46, 707 (1966). 

11. J. Gasser and H.Leutwyler, Phys. Rep. 87, 77 (1982). 

12. A. W. Thomas and W. Weise, The Structure of the Nu- 
cleon, (WILEY- VCH, Berlin, 2001) chapter 7. 

13. S.Weinberg, Physica 96A, 327 (1979). 

14. J. Gasser and H.Leutwyler, Ann. Phys. 158, 142 (1984). 

15. J. Gasser and H.Leutwyler, Nucl. Phys. B 250, 465 
(1985). 

16. G.Ecker, Prog. Part. Nucl. Phys. 35, 1 (1995). 

17. V.Bernard, N. Kaiser, and Ulf-G. Meifiner, Int. J. Mod. 
Phys. E 4, 193 (1995). 

18. S. Scherer, Introduction to Chiral Perturbation The- 
ory, Adv. Nucl. Phys. 27, 277 (2003) | |arXiv:hep-| 
ph/0210 398vl). 

19. V.Bernard, Prog. Part. Nucl. Phys. 60, 82 (2008). 

20. T. Becher and H.Leutwyler, J. High En. Phys. 06, 017 
(2001). 

21. J. Gasser, H.Leutwyler, and M. E. Sainio, Phys. Lett. B 
253, 252 (1991). 

22. M. E. Sainio, Pion-nucleon a-term - a review, in Proc. of 
the 9 th Symp. on Meson-Nucleon Physics and the Struc- 
ture of the Nucleon (MENU'01 ), Washington D. C, 2001, 
tvN newsletter 16, 138 (2002), ISSN 0942-4148, and ref- 
erences therein. 

23. M. L. Goldberger, H. Miyazawa and R. Oehme, Phys. Rev. 
99, 986 (1955). 

24. T. E.O.Ericson, B.Loiseau and A.W.Thomas, Phys. 
Rev. C 66, 014005 (2002). 

25. V.V. Abaev, P. Metsa, and M. E. Sainio, Eur. Phys. J. A 
32, 321 (2007). 

26. V.Bernard, N. Kaiser, and Ulf-G. Meifiner, Phys. Lett. B 
383, 116 (1996). 

27. M. A.Kovash et al., irN newsletter 12, 51 (1997). 

28. O. Hanstein, D. Drechsel, and L. Tiator, tvN newsletter 
12, 56 (1997). 

29. D.F. Measday, Phys. Rep. 354, 243 (2001) . 

30. T. Gorringe and H. W. Fearing, Rev. Mod. Phys. 76, 31 
(2004). 

31. V.A. Andreev et al, Phys. Rev. Lett. 99, 032002 (2007). 

32. D.Sigg et al., Nucl. Phys. A 609, 310 (1996). 

33. V. E. Lyubovitskij and A. Rusetsky, Phys. Lett. B 494, 9 
(2000). 

34. J. Gasser et al, Eur. Phys. J. C 26, 13 (2002). 

35. P. Zemp, Proc. of Chiral Dynamics 2003, p. 128, 
Bonn, Germany, September 8-13, 2003, arXiv:hep- 
ph/0311212^ 1; Pionic Hydrogen in QCD+QED: Decay 
width at NNLO. Ph. D. Thesis, University of Bern (2004). 

36. G.Rasche and W. S. Woolcock, Nucl. Phys., A 381, 405 
(1982). 

37. J. Spuller et al, Phys. Lett. 67 B, 479 (1977). 

38. G. C.Oades, G.Rasche, W.S. Woolcock, E. Matsinos, 
A.Gashi, Nucl. Phys. A 794, 73 (2007). 

39. T. E. O. Ericson and W. Weise, Pions and Nuclei, (Claren- 
don, Oxford, 1988) chapter 9. 

40. ibidem, chapter 4. 

41. A.W.Thomas and R. H. Landau, Phys. Rep. B 58, 121 
(1980). 

42. A. Deloff, Fundamentals in Hadronic Atom Theory, 
(World Scientific, London, 2003) chapter 15. 



43. D.Sigg et al, Nucl. Phys. A 609, 269 (1996). 

44. D. Chatellard et al, Phys. Rev. Lett. 74, 4157 (1995); 
Nucl. Phys. A 625, 855 (1997). 

45. PHauser et al, Phys. Rev. C 58, R1869 (1998). 

46. H.-Ch. Schroder et al, Eur. Phys. J C 21, 473 (2001). 

47. S.Weinberg, Phys. Lett. B 295, 114 (1992). 

48. V. Baru and A. E. Kudryavtsev, Phys. of At. Nucl. 60, 
1475 (1997). 

49. S.R. Beane, V.Bernard, T.-S. Lee, and U.-G. Meifiner, 
Phys. Rev. C 57, 424 (1998). 

50. V. E. Tarasov, V. V. Baru, and A. E. Kudryavtsev, Phys. 
At. Nucl. 63, 801 (2000). 

51. A. Deloff, Phys. Rev. C 64 (2001) 065205. 

52. N. Kaiser, Phys. Rev. C 65, 057001 (2002). 

53. B. Burasoy and H. W. Grieshammer, Int. J. Mod. Phys. E 
12, 65 (2003). 

54. S. R. Beane, V. Bernard, E. Epelbaum, U.-G. Meifiner, 
and D.R.Phillips, Nucl. Phys. A 720, 399 (2003). 

55. M. Doring, E.Oset, and M.J.Vicente Vacas, Phys. Rev. 
C 70, 045203 (2004). 

56. B. F. Irgaziev and B. A. Fayzullaev, arXiv:hep- 
ph/0404203vl (2004). 

57. U.- G. Meifiner, U. Raha, and, A. Rusetsky, Eur. Phys. J. 
C 41, 213 (2005); Eur. Phys. J. C 45, 545 (2006). 

58. U.- G. Meifiner, U. Raha, and, A. Rusetsky, Phys. Lett. B 
639, 478 (2006). 

59. M. Hoferichter, B. Kubis and U.-G. Meifiner, Phys. Lett. 
B 678, 65 (2009); Nucl. Phys. A 833,18 (2010). 

60. V. Baru et al, Phys. Lett. B 694, 473 (2011) (arXivmucl- 
th/1003.4444v2). 

61. M. Hoferichter, B. Kubis and U.-G. Meifiner, PoS CD09, 
014 (2009) (arXiv:hep-ph/0910.0736). 

62. J.Hufner, Phys. Rep. 21, 1 (1975). 

63. V.Lensky et al, Eur. Phys. J. A 27, 37 (2006). 

64. V.Lensky et al, Phys. Lett. B 648, 46 (2007). 

65. V.L.Highland et al, Nucl. Phys. A 365, 333 (1981). 

66. D.W.Joseph, Phys. Rev. 119, 805 (1960). 

67. R.MacDonald et al, Phys. Rev. Lett. 38, 746 (1977). 

68. PSI proposal R-98-01, www2.fz-juelich.de/ikp/exotic- 
atoms. 

69. D. Gotta et al, Lect. Notes Phys. 745, 165 (2008). 

70. Th. Strauch et al, Phys. Rev. Lett. 104, 142503 (2010). 

71. F.J. Hartmann, Proceedings of Physics of Exotic Atoms 
on Electromagnetic Cascade and Chemistry, Erice, Italy, 
1989 (Plenum Press, New York, 1990) p. 23 and p. 127, 
and references therein. 

72. J.S.Cohen, Rep. Prog. Phys. 67, 1769 (2004). 

73. E.Borie and M.Leon, Phys. Rev. A 21, 1460 (1980). 

74. M.Leon and H. A. Bethe, Phys. Rev. 127, 636 (1962). 

75. R. al Hassani et al, Z. Phys. A 351, 113 (1995). 

76. T. S. Jensen and V. E. Markushin, Eur. Phys. J. D 19, 165 
(2002). 

77. T. S. Jensen and V. E. Markushin, Eur. Phys. J. D 21, 261 
(2002). 

78. T. S. Jensen and V. E. Markushin, Eur. Phys. J. D 21, 271 
(2002). 

79. L.Bracci and G. Fiorentini, Nuovo Cim. A 43, 9 (1978). 

80. R.Pohl et al, Phys. Rev. Lett. 97, 193402 (2006). 

81. J.B. Czirr et al, Phys. Rev. 130, 341 (1963). 

82. A. Badertscher et al, Europhys. Lett. 54, 313 (2001), and 
references therein. 

83. D. S.Covita et al, Phys. Rev. Lett. 102, 023401 (2009). 



Th. Strauch et al: Pionic Deuterium 



19 



84. M. Hennebach, Precision measurement of ground state 
transitions in pionic hydrogen, Ph.D. thesis, Universitat 
zu Kohi (2003), [htFp://kups.ub.urii-koeln.de/744/[ 

85. V. N. Pomerantsev and V. P. Popov, Phys. Rev. A 73, 
040501(R) (2006). 

86. T.S.Jensen, V.N.Pomerantsev, and V.P.Popov, 
arXiv:nucl-th/0712.3010vl (2007). 

87. V. P. Popov and V. N. Pomerantsev, arXiv:nucl- 
th/0712.3111vl (2007). 

88. D.Taqqu, AIP Conf. Proc. 181, 217 (1989). 

89. S. Jonsell, J. Wallenius, P. Froelich, Phys. Rev. A 59, 3440 
(1999). 

90. R.Pohl, Hyperfine Int. 193, 115 (2009). 

91. E.Lindroth, J. Wallenius, S. Jonsell, Phys. Rev. A 68, 
032502 (2003); Phys. Rev. A 69, 059903(E) (2004). 

92. S.Kilic, J.-PKarr, L.Hilico, Phys. Rev. A TO, 042506 
(2004). 

93. T. L. Trueman, Nucl. Phys. 26, 57 (1961). 

94. E.Lambert, Helv. Phys. Acta 42 (1969) 667; ibidem 43, 
713 (1970). 

95. J. Mitroy and I. A. Ivallov, J. Phys. G 27, 1421 (2001). 

96. D.Eiras and J.Soto, Phys. Lett. B 491, 101 (2000). 

97. A. H. Rosenfeld, Phys. Rev. 96, 139 (1954). 

98. A.Reitan, Nucl. Phys. B 11, 170 (1969). 

99. H. Machner and J.Niskanen, Nucl.Phys. A 776, 172 

(2006) . 

100. K. Bruckner, R. Serber, and K.Watson, Phys. Rev. 81, 
575 (1951). 

101. V. Baru, C.Hanhart, and A. Rusetski, priv. comm. 

102. A.Filin et al, Phys. Lett. B 681, 423 (2009). 

103. E. Byckling and K. Kajantie, Particle Kinematics, John 
Wiley and Sons, London (1973) chapter II. 6. 

104. L.M.Simons, Physica Scripta T22, 90 [1988). 

105. L.M.Simons, Hyperfine Int. 81, 253 (1993). 

106. D. S. Covita et al, Rev. Scient. Instr. 79, 033102 (2008). 

107. J. Eggs and K. Ulmer, Z. angew. Phys., 20(2), 118 (1965). 

108. G. Zschornack, Nucl. Instr. Meth. 200, 481 (1982). 

109. N.Nelms et al, Nucl. Instr. Meth. A 484, 419 (2002). 

110. G.Basile et al, Phys. Rev. Lett. 72, 3133 (1994). 

111. R. Deslattes et al, Rev. Mod. Phys., vol. 75, no. 1, 35 
(2003). 

112. P.Indelicato et al, Rev. Sci. Instrum. 77, 043107 (2006). 

113. D. F. Anagnostopulos et al, Nucl. Instr. Meth. A 545, 
217 (2005). 

114. M. Trassinelli et al, J. Phys., Conf. Ser. 58, 129 (2007). 

115. D. S. Covita, High-precision spectroscopy of the 3p - Is 
transition in muonic hydrogen, Ph. D. thesis, University 
of Coimbra (2008), [http~77 hdl.handl e.net/10316/7521| 

116. M. Sanchez del Rio and R. J. Dejus, Proc. SPIE Int. Soc. 
Opt. Eng. 3448, 246 (1998); ibidem 5536, 171 (2004); 
|http: / /www, esrf ,eu7 computing / scientific /xo p2 . 1 1 

117. Th. Strauch, High-precision measurement of strong- 
interaction effects in pionic deuterium, Ph. D. the- 
sis, Universitat zu Koln (2009), |http://kups.ub.uni-| 
koeln.de/2813/ 

118. WM. J. Veigele, At. Data Tables 5, 51 (1973). 

119. M. O.Krause and J.H.Oliver, J. Phys. Chem. Ref. Data 
8, 329 (1979). 

120. T. Mooney, Argonne National Laboratory, priv. comm. 

(2007) . 

121. B.L.Henke, E. M. Gullikson, and J.C.Davies, At. Data 
Nucl. Data Tables 54, 181 (1993). 

122. C. T. Chantler, J. Phys. Chem. Ref. Data 24, 71 (1995). 



123. F. N. Chukhovskii, G. Holzer, O. Wehrhan, and 
E.Forster, J. Appl. Cryst. 29, 438 (1998). 

124. F. Cembali et al, J. Appl. Cryst. 24, 424 (1992). 

125. F. James and M. Roos, Comput. Phys. Commun. 10, 343 
(1975). 

126. U. C. Bergmann and K. Riisager, Nucl. Instrum. Meth. A 
489, 444 (2002). 

127. C. Amsler et al (PDG 2008), Phys. Lett. B 667, 1 (2008). 

128. R.Pohl et al, Nature 466, 213 (2010). 

129. C. G. Parthey et al, Phys. Rev. Lett. 104, 233001 (2010). 

130. Y.Lu and R. Rosenfelder, Phys. Lett. B 319, 7 (1993); 
Phys. Lett. B 333, 564(E) (1994). 

131. W. Leidemann and R. Rosenfelder, Phys. Rev. B 51, 427 
(1995). 

132. J. Ahrens et al, Eur. Phys. J. A 23, 113 (2005). 

133. R.Koch, Nucl. Phys. A 448, 707 (1986). 

134. J.Bailey et al, Phys. Lett. 50 B, 403 (1974). 

135. E.Bovet et al, Nucl. Instr. Meth. 190, 613 (1981). 

136. E.Bovet et al, Phys. Lett. 153 B, 231 (1985). 

137. C.Hanhart, priv. comm. (2010). 

138. F. S. Crawford and M. L. Stevenson, Phys. Rev. 97, 1305 
(1955). 

139. C.M.Rose, Phys. Rev. 154, 1305 (1967). 

140. C.Richard-Serre et al, Nucl. Phys. B 20, 413 (1970). 

141. D.Aebischer et al, Nucl. Phys. B 108, 214 (1976). 

142. D.A.Hutcheon et al, Phys. Rev. Lett. 64, 176 (1990); 
Nucl. Rev. A 535, 618 (1991). 

143. B.G.Ritchie et al, Phys. Rev. Lett. 66, 568 (1991). 

144. PHeimberg et al, Phys. Rev. Lett. 77, 1012 (1996). 

145. M.Drochner et al, Nucl. Phys. A 643, 55 (1998). 

146. D. S. Koltun and A. Reitan, Phys. Rev. 141, 1413 (1966). 

147. I.R.Afnan and A.W.Thomas, Phys. Rev. C 10, 109 
(1974). 

148. C.Fayard, G.H. Lamot, and T. Mizutani, Phys. Rev. 
Lett. 45, 524 (1980). 

149. C.J.Horowitz, Phys. Rev. C 48, 2920 (1993). 

150. J. A. Niskanen, Phys. Rev. C 53, 526 (1996). 

151. A.N.Ivanov et al, Z. Phys. A 358, 81 (1997). 

152. S.G.Eckstein, Phys. Rev. 129, 413 (1963). 

153. J.-F. Germond and C.Wilkin, J. Phys. G 14, 181 (1988). 

154. D. Gotta et al, Phys. Rev. C 51, 469 (1995). 

155. E.Daum et al, Nucl. Phys. A 589, 553 (1995). 

156. J. Roginsky and C. Werntz, Phys. Rev. C 40, 2732 (1989). 

157. D.G.Long, M.Sternheim, and R. R. Silbar, Phys. Rev C 
26, 586 (1982). 

158. G. A. Miller and P. Sauer, Phys. Rev. C 44, R1725 (1991). 

159. H.O.Meyer et al, Nucl. Phys. A 539, 633 (1992). 

160. C.Hanhart, Phys. Rep. 397, 155 (2004). 

161. V.Baru et al, Phys. Rev. C 80, 044003 (2009). 

162. C.Hanhart and N. Kaiser, Phys. Rev. C 66, 054005 
(2002). 



